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Abstract 

We determine the average number d (N, K), of A^C-Kauffman networks 
that give rise to the same binary function. We show that, for N ^> 1, there 
exists a connectivity critical value K c such that i?(iV, K) pa e^ (ip > 0) for 
K < K c and i?(iV, K) pa 1 for K > K c . We find that K c is not a constant, 
but scales very slowly with N, as K c pa log 2 log 2 (2A/ln2). The problem of 
genetic robustness emerges as a statistical property of the ensemble of NK- 
Kauffman networks and impose tight constraints in the average number of 
epistatic interactions that the genotype-phenotype map can have. 
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1. Introduction 



A^-Kauffman networks, also known as TVi^-Kauffman cellular automata, 
were proposed in 1969 as models for the study of gene regulation 1,2 . Since 
then, their rich dynamical behavior has motivated many studies, and as very 
general models with few parameters their applications have been extended 
to several complex systems. Of particular interest in this work is their appli- 
cation for modeling the genotype-phenotype map 3 . For an excellent review 
on several applications see also Ref. 4. 

An A^if-Kauffman network consists of N Boolean variables (or bits) 
Si(t) G {0, 1}, with % — 1, . . . , N, that evolve deterministically in discretized 
time t = 0, 1, 2, . . . according to Boolean functions on K of these variables 
(the inputs) at the previous time t — 1. For every Si, a Boolean function 
fi is chosen at random and independently from a given distribution in all 
the possible Boolean functions with K inputs. Also, for every Si, K inputs 
are randomly selected from a uniform distribution among the iV Boolean 
variables of the network. The selection process may be done with allowed 
repetitions; this is to say, some inputs might be identical and K could be 
bigger than iV 5 . Or without repetitions; namely, all inputs are different and 
K < N . In this work we adopt the latter alternative, because it is more nat- 
ural from the biological point of view, and it suits better for our calculation 
purposes. It is important to bear in mind this distinction (repetition vs. no 
repetition) when contrasting our results with those from other approaches. 

Once the K inputs and the function fi for each of the N variables have 
been selected, a particular AfAT-Kauffman network has been constructed, and 
evolves deterministically in discrete time t according to the rules 

Si(t + l) = fi (S n (t), S i2 (t), S iK (t)) , i — l,...,N, (1) 

where i a ^ ip, for all a, (3 — 1, 2, ... , K, with a ^ f3, because all the inputs, 
while random, are different. 

Lets denote C% the set of different A^-Kauffman networks that might 
be built up with the aforementioned process, for given N and K. Since each 
element of C^- is a deterministic dynamical system evolving in a finite phase 
space, with 2 N states, its dynamics eventually settles into a cycle. One can 
think the system as composed of attraction valleys, each with one cycle (or 
attractor), whose number might go from just 1 to 2 N , and with cycle lengths 
varying from 1 for a punctual attractor, to 2^ when a single cycle traverses 
the whole phase space. The behavior of a typical A^-Kauffman network is 
in between these two extreme cases. 

Determining the distribution in of the number of attractors and their 
size for general and K, is a difficult and challenging problem. Several 
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numerical simulations and analytical methods have been used to approach 
it 5 " 10 . However, significant advances have been obtained only for the so- 
called random map model; i.e., when the Boolean functions are taken from 
a uniform distribution and K = N (K = oo in Ref. 5, since they permit 
repetition in the inputs). The first results obtained from this model go as far 
as the 50 's, when accurate formulas were found for the distribution of attrac- 
tors in the context of random functional graphs 11 . Decades later, Derrida 
and Flyvbjerg studied it in the context of statistical mechanics with an inter- 
est on its applications to cellular automata 5 . We recently have founded an 
asymptotic formula for the statistical distribution of the number of connected 
components in the random map model, deriving it from a new combinato- 
rial expression 12 . In Ref. 13 we also, furnished both: exact and asymptotic 
formulas, for several measures that help to understand the connectivity of 
random functional graphs; such as cycle and trajectory lengths, expected 
number and size of attraction valleys, and the like. 

The general case K < N were treated by mean field analysis taking 
N — > oo 10 . The results indicate that when the Boolean functions are ex- 
tracted from a distribution such that fi is 1 or with probability p or 1 — p, 
respectively, there is a bifurcation of the dynamics of the elements of at 
the critical connectivity value 

K* = — 7^ r , (2) 

with p — 1/2 corresponding to the uniform distribution considered in this 
paper. For K < K*; there is an ordered phase, where small perturbations 
die out. When K > K*, the phase is chaotic and small perturbations spread 
exponentially through the network. While, in the critical case K = K*, the 
evolution is mainly governed by fluctuations and has been qualified of being 
neither ordered, nor chaotic 4 . 

The extreme cases K — 1 and K = 2 have been approached analytically 
for p — 1/2. The former was first studied by Flyvbjerg and Kjaer 6 , and it 
has been recently founded that, when the constant functions tautology and 
contradiction are excluded from consideration, the number of attractors and 
their length grow super polynomially (faster than any power law) in N 7 . In 
the case of K = 2 a super polynomial behavior for the number of attractors 
was also founded 8 . 

In this work we calculate; for given N and K, the exact value of the 
average number d (N, K) of different iVX-Kauffman networks that give rise 
to the same binary function. For that scope, we define the function \1/ : £^ — > 
G2 N 'i where Q 2 N is the set of binary functions in iV binary variables. Then, 
we calculate the values of the cardinalities and ! °f £-k an d 

the set of binary functions, ^ (£-k) — ^2 N that they generate, respectively. 
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Our findings, show that, for N 1 the asymptotic formula for i?(iV, K), 
behaves so that there exists a critical value K c of the connectivity, such 
that: ${N,K) « e ^ K)N for K < K c , with ip(K) > 0; so that many NK- 
Kauffman networks generate the same binary function. And, i?(iV, if) « 1 
for K > K c , indicating that almost any binary function is generated by a 
different iVi^-Kauffman network. Furthermore, the value of K c turns out 
to be not a constant, but to grow very slowly as the double logarithm of 
N. Important to remark, is the fact that, K c does not signal a transition 



from a regular to a chaotic behavior for the elements of C^, as K* given by 
(2) does. Instead, it shows an abrupt change in the injectivity of the map 



iVif -Kauffman networks play an important roll in applications to ge- 
netics for modeling the genotype-phenotype map, represented by \P 3 . The 
genotype, carries all the necessary information to create a living organism; 
the phenotype. The genotype is mainly composed of DNA that is a dou- 
ble chain of base-pares (Adenine-Thymine and Cytosine-Guanine) which in 
turn constitutes the genes 14 . In this context, an A^-Kauffman network, 
represents the genotype, while their attractors in \I/ {£>k) ] the phenotype of 
the alternative cells types 1,15 . A conspicuous observation in the theory of 
natural selection is the robustness of the phenotype against mutations in the 
genotype 16_18 . Natural radiation in the environment changes the genotype 
by making mutations on DNA-bases. Random mutations and recombination 
of the genotype by mating, constitute the essential engine of species evolu- 
tion. The effect of random mutations range from having no effect at all, to 
complete damage in the phenotype. The reasons of this phenomena are still 
not completely understood. One hypothesis is that it provides selective ad- 
vantages to the phenotype since an organism with a damaging mutation will 
be at disadvantage in evolution 18 . So, natural selection must have favored 
organisms with a mechanism that prevents mutations from accumulating in 
the gene, ensuring that useful genes remain in the genome 14 . 

A mechanism proposed to explain genetic robustness is due to the finding 
of the existence of genetic redundant material 14 . This implies that not all 
genes are essential. Experiments to estimate how many genes are actually 
essential were carried out by induced mutations 19 . The results varied among 
the different organisms under study, but in all cases; the estimates showed 
that more than 50% of the genes are not essential. However, recent studies 
on gene-deletion at genome-scale, revealed that thousands of genes whose 
deletions had no detectable effect in the phenotype were single-copy i. e. they 
had no duplicate in the genome 17 . Other mechanism is the discovery by 
modern genetics of the epistatic effects. This refers to effects that a gene 
may have on the phenotype, that strongly depend on the levels of expression 
of other unrelated genes in the genotype 16 . For example, studies in yeast 
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revealed that up to 50% of mutations, almost do not affect their fitness due to 
the epistatic compensations 20 . Today there is some agreement that genetic 
robustness, emerges as a mixture, between genetic redundancy and epistatic 
buffering 18 . 

As we shall see, our mathematical findings show that genetic robust- 
ness emerges in the genotype-phenotype map modeled by A^-Kauffman 
networks, as a consequence of their statistical properties: with the K connec- 
tions playing the roll of the average number of epistatic interactions among 
the genes. Our calculations impose tight restrictions on the values that K 
may have when N attains values in the range of the known number of genes 
that living organisms have. 

The paper is organized as follows: In Sec. 2 we establish a mathemati- 
cal correspondence between binary functions and functional graphs through 
a bijection This allows us to define the function ^ and show that 

•&(N,K) = #£^/# , l/ (£k)- The exact computations of the cardinalities 
#£x> and (£^); are carried out by combinatorial methods in Sec. 3. 
In Sec. 4 the asymptotic behavior of •§ (N, K) is studied and expressions for 
K C (N) and AK C (N) (the width of the transition) are found. Finally, in Sec. 5 
we set our conclusions showing that: the asymptotic behavior of •d (N, K) im- 
poses the restriction K < 3, in the average number of epistatic interactions, 
in order that, genetic robustness emerges in the ensemble of A^-Kauffman 
networks. 



2. Mapping Af/f-Kauffman Networks into Binary Functions 



Binary functions can be represented by means of functional graphs 12 > 13 . For 
a given positive integer n, an n-functional graph is composed: by the set 
P n — {1, . . . , n}, whose elements are called nodes or vertices, and a function 

9 : Pn - Pn- 

Functional graphs can be represented graphically, and it is this graphical 
representation that helps to understand many of their properties. Each node 
is indicated by a point, and an arrow from node i to node j is drawn whenever 
g(i) = j. As an example, Figure 1 shows a 12-functional graph with three 
connected components (i.e., three attractors). For clarity, nodes are depicted 
here as small disks. 



5 




Figure 1: A 12-functional graph with three connected components. 

Note that, in an n-functional graph, although exactly one arrow goes out 
from each of its nodes, several arrows might be directed towards any of them. 
Denote Q n the set of all possible n-functional graphs; its cardinality is well 
known and given by j^Q n = n n 13 . The problem of determining in Q n the 
distribution of cycle lengths and number of connected components (number of 
attractors) has been undertaken by several authors since the early 50's n ~ 13 . 

Let Q N = { S = (Si, . . . , S N ) \ Si = 0,1, for % = 1, . . . , TV } be the set 
of 2 N binary vectors with N components which is mapped by A/7f-Kauffman 
networks (/ : Qn — > Qn) through (1). There is a well-known bijection 
07V : Qn — ► P2 N given by: 

N 

cf> N (S) = l + J2 2i ~ 1 Si- (3) 

i=i 

Clearly, S is nothing else than the binary decomposition of 0a? (S) G P 2 n. 
Now, the following diagram commutes: 



0A 

P 2 N 



f 



g 



0A 

P 2 N 



and assigns to each binary function / the functional graph g = <f) N o f o <f) N . 
This defines the function 

^:C N K ^^{C N K ) CQ 2N . (4) 
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In the case K = N we get from (1) that £^ coincides with the set of all 
possible binary functions from Qn — > Qn- So, from the commuting diagram 
it happens that £^ = Q 2 N ■ 

The average number of A^-Kauffman networks that \1/ maps to the same 
binary function is directly expressed in terms of the cardinalities of the inverse 
image sets tt" 1 (g) = {/ G C% | (/) = g) as 

where 

r 1 ( ff )nr 1 ( 9 ') = $ Vg^g'- (6) 

Furthermore £^ may be decomposed as 

c%= U 

so that, due to (6) 

= £ #^ («?) . 

Substituting back into (5) we obtain 



* W = ^ (7) 



for the average number of A/if-Kauffman networks that give rise to the same 
binary function. 



3. The Cardinalities of and (C%) 



For each of the A^ Boolean variables there are 2 2 different Boolean func- 
tions with K connections; moreover, as there are (^) different ways to make 
the connections without replacement, the total number of A^-Kauffman 
networks is 

N 

(8) 
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Now, let Ax = {ijf(m, i)} denote the 2 K x K binary matrix whose m-th 
row encodes the binary decomposition of m G P 2 k , that is, 

K 

m = l + ^A K (m,i)2 i ~ 1 . 
i=i 

As an example, A 3 is shown in Figure 2, at left. 



.(ii) 



.(12) 



.(13) 





1 

1 

1 1 

1 
1 1 
1 1 

1 1 1 



" 





" 




" 





1 " 




" 





" 





1 


1 







1 


1 







1 





1 










1 










1 





1 


1 


1 


1 




1 


1 


1 




1 


1 


1 



Figure 2. 



There are (^) possible fT-connection sets 

c k } = Oi) *2, • • • , ix} Q Pn, with a = 1, 



(9) 



and hence without loss of generality we can take i\ < i 2 < • • • < Ik- To each 



we associate a i^-connection map C K a> : Qn —> such that 



*(a) 



(S) — C K a> (Si,..., Sn) — (S^, . . . , S iK ) . 
A i^-Boolean function 



(10) 



(11) 



is completely determined by a if-truth table T K , that consists of a 2^ x 
(i^ + 1) binary matrix of the form 



T K = [A K h K ] , 



(12) 
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where the vector 

b K = [<7i,<7 2 , • • -,<72k] (13) 

is expressed in column form in (12), and its entries correspond to the 2 K 
images of the function (11). 

As any X-truth table (12) has 2 K rows there are as much as 2 2K possible 
if-truth tables, i.e., the total quantity of binary vectors (13). The Boolean 
functions (11) have been classified by Wolfram's notation according to their 
decimal number fi given by 21 

/i=l + ^2 s -V s . (14) 

s=l 



Let us number the vectors (13) as well as the If -truth tables (12) accord- 
ing to (14) by adding the superscript (//), i.e., and 



K — 



A K 



with /i = 1, . . . ,2 K . Thus, for instance, the 2-truth tables T 2 , T 2 , and 
T 2 (13) in Figure 2. 

The connection function (10) projects Qn on Qk, enabling us to define 

jected iV-Boolean functions b^ 1 "" 
defined through their iV-truth tables 



projected iV-Boolean functions = b^ o C*^ a \ that are completely 



T 



(*0(«) 



A' 



A 



N 



D K ° °K 



(15) 



Within this notation, an MT-Kauffman network consists of a set 




and an evolution rule (equivalent to (1)): 

s i (t + i) = b ( p ) ocz*(s{t)), 

where some of the indexes ctj and ^ may be the same for different values of 
i; and a 2 Ar -functional graph g e G 2 n is associated through \1/ in (4). 

An example is in point. Let N = 3, K = 2, and consider the 2-connection 
sets = {1,2}, Cf ] = {1,3}, cf } = {2,3}, and the 2-truth tables T 2 (11) , 
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T 2 (12) , T 2 {13) (see Figure 2). We can construct an iV-fT-Kauffman network 
consisting of the projected tables {T, 



; il)(1) ,T 3 (12)(2) ,r 3 (11)(3) }, that gives rise 
to an 8-functional graph with two connected components as in Figure 3. 

Note that taking the projected tables {T 3 (13)(3) , T 3 (12){2) , T 3 (11){2) } the same 
8-functional graph appears. 



010 




101 
• — 



111 



"Q 



100 110 001 
Figure 3. The nodes are labelled here in binary form, according to (3). 



Now, for a = 1, . . . , (^), let B^{a) denote the 2^ x 2 2K binary matrix 
whose fi-th column results from the application of the projected TV-truth 



table T^ )(a) given by (15), that is: 



B N K {a) 



*(<*) 



>K 



oC 



*{a) 
K 



To clarify how matrices B^(a) are built up, consider for example N = 3 
and K = 2. There are 2 4 = 16 truth tables, numbered according to (14), 
and shown below in compact form. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


00 





1 





1 





1 





1 





1 





1 





1 





1 


01 
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1 
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1 
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1 








1 


1 


10 
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1 


1 


1 
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1 


1 


1 


11 
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1 


1 


1 


1 


1 


1 



Moreover, there are ( 3 ) = 3 possible ways to choose two columns of A 3 . 
Thus, assuming that columns 1 and 2 of A3 form the connection set number 
1 we get Bl(l): 



10 



^3 
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10 


11 


12 


13 


14 


15 


16 


000 





1 





1 





1 





1 





1 





1 





1 





1 


001 





1 





1 





1 





1 





1 





1 





1 





1 


010 








1 


1 








1 


1 








1 


1 








1 


1 


Oil 








1 


1 








1 


1 








1 


1 








1 


1 


100 














1 


1 


1 


1 














1 


1 


1 


1 


101 














1 


1 


1 


1 














1 


1 


1 


1 


110 


























1 


1 


1 


1 


1 


1 


1 


1 


111 


























1 


1 


1 


1 


1 


1 


1 


1 



Taking columns 1 and 3 of A 3 as the connection set number 2 yields matrix 
Bl(2): 



A 3 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


000 





1 





1 





1 





1 





1 





1 





1 





1 


001 








1 


1 








1 


1 








1 


1 








1 


1 


010 





1 





1 





1 





1 





1 





1 





1 





1 


Oil 








1 


1 








1 


1 








1 


1 








1 


1 
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1 
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1 
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1 
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1 


1 


1 



Finally, matrix £?f(3) results from columns 2 and 3 of A 3 : 



A 3 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


000 





1 





1 





1 





1 





1 





1 





1 





1 


001 








1 


1 








1 


1 








1 


1 








1 


1 


010 














1 


1 


1 


1 














1 


1 


1 


1 


011 


























1 


1 


1 


1 


1 


1 


1 


1 


100 





1 





1 





1 





1 





1 





1 





1 





1 


101 








1 


1 








1 


1 








1 


1 








1 


1 


110 














1 


1 


1 


1 














1 


1 


1 


1 


111 


























1 


1 


1 


1 


1 


1 


1 


1 



Let 

The number {£-k) °f different functional graphs that can be gener- 
ated by the iVi^-Kauffman networks, equals the N-ih power of the number 
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c of different columns belonging to the matrices of S^. Moreover, 




where r is the number of 'redundant' columns (that is, the ones that are 
repeated) in the elements of S^. As no two columns of A N are the same, 
redundant columns arise solely from the following functions: 



• The tautology: bx {C* K (S)) = 1, for any S G fi/v. As there is one 
tautology (column 16 in the last example) in each of the (^) matrices, 
we get as much as (^) — 1 redundant columns. 

• The contradiction: bx (C* K (S)) = 0, for any S G VLn- With one con- 
tradiction (column 1 in the last example) in each matrix, we come also 

to (^) — 1 redundant columns. 

• The identity: b K (S il , . . . , S iK ) = S iv for / = \,...,K. Observe that 
no two columns of An are the same. Hence, for each j = 1, . . . ,N, 
one replica of column j of An can be found in precisely (^Z\) matri- 
ces, yielding (^~*) — 1 redundant columns. Hence there is a total of 
N [(^Ii) — l] redundant columns accounting for the identity function, 

(in the last example, columns 11 of -Bf(l) and 13 of £>f (3) replicate 
column 2 of A 3 ). 

• The negation: b K (S ii: . . . , S iK ) = 1 — S iv for I = 1,...,K. This 
function leads also to N [(^~*) — l] redundant columns, because one 
complement of column j of An can be found in each of (^Z\) matrices. 



Summing up we get 



-2 + 2N 



N — 1 
K - 1 



-2N, 



as the total number of redundant columns in S^. Substituting in (16) and us- 

ing the identity (£:}) = #(£), yields c = {tf K -2(K+ 1)) © +2 (N + 1). 

Thus, the number of different functional graphs generated by NK-K&u&m&n 
networks amounts to 



(2 2X -2(K+1)) +2(N + l 



N 



(17) 
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4. The Asymptotic Expansion of -d (N, K) 



From equations (7), (8) and (17) we get for the reciprocal of d (N, K) that: 
0" 1 (N,K) — {1 — (p (K) [1 - i N (K)]} N , (18) 



where 



and 



N + 1 /iV\ _1 



The function 1;n(K) is unimodal on the interval 1 < K < N, where it 
attains its maximum at K = N, yielding £n{N) = 1, and its minimum at 

K = N/2, where £ N (N/2) ~ C . 

For fixed iV, though TV >> 1, we obtain the following: 

By means of Stirling's approximation in the factorials involved in (^) \ 
it happens that £ N (K) ~ C^l/iV^- 1 ) for if ~ (9(1), and &v(iV - m) ~ 
C(l/iV m ) for m ~ 0(1). Thus, £jv(if) ~ o(l) over 1 < AT < iV. Hence, since 
</?(if) <C 1 for if > 1, t? -1 (7V, if) exhibits the following asymptotic behavior 

&-\N,K) w {l-^J^^we-^. (20) 

Figure 4 shows an excellent agreement for the functions (18) and (20) for 
just iV = 10. There are two main asymptotic regimes: 

i) tT 1 (N,K) « 1 for 7VV(iO < 1, (21.a) 

and 

it) i? -1 (JV, if) ~ for iV>(if) > 1, (21.6) 

The critical point if c defining the transition region; where the regime changes 
is given by setting -d~ l (N,K c ) = 1/2 in (20), yielding the transcendental 
equation 

2^n2-ln(if c +l)=ln(j^, (22) 



with solution 



r, , , C*N\ , n (\n\a\nN\ 

^«io g2 io g2 - ) + o Kr^l- < 23 > 
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Figure 4: The graph of $ 1 (N, K) as a function of K and its asymptotic approximation 
(dashed); for N = 10. 



The leading term corresponds to the solution obtained by neglecting the 
second term in the left hand side of (22). Equation (23) gives K c as a very 
slowly growing function on N. To estimate the width of the transition region 
let us expand ^^(N, K) in Taylor series up to the first order in K — K c . 
From (20) we obtain: 

tf" 1 (A^K) ^ i [1 - AV (K c ) (^ - , for \K-K C \<^1. (24) 

The width AK C of the transition region is then given by: 
AK C = K 1 -K = 



Ay (K c y 

where K and K l are such that tf -1 (N,K ) = and (iV, = 1, in 
(24). From (19) we get 

AK C = - 2 J&±V O(JL), (25) 

[2^(K c + l)(ln2) 2 - 1] In 2 \lnN J K ! 
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that is smaller than the absolute error in (23). 

Summing up the results from (20) and (21) we have, for the asymptotic 
regimes outside of the transition region (of width AKc), that: 

( e^ N > 1 for K < K c — AK C 
&(N,K)f*l . (26) 

{ 1 + (p(K) N « 1 for K > K c + AK C 



5. Conclusion 



We have calculated an exact formula (18) for the average number d (N, K) 
of ATif-Kauffman networks, that are mapped by \I/ [defined by (4)] onto the 
same binary function. The asymptotic expression (26) for d (N, K) shows an 
abrupt change of regime at the critical value K c (23), that grows with N as 
a double logarithm. The width of the transition AK C (25), becomes small as 
0(l/\nN). 

In genetics : -> * (£%) C Q 2 n may be used for modeling the 

genotype-phenotype map 3 . The genotype is represented by a A^-Kauffman 
network with N Boolean variables; Si (t) representing the expression of the 
i-th gene at time t within some developmental process. While a gene's ex- 
pression, could be much more complex than just to be described by binary 
values; ATi-f-Kauffman networks have enough mathematical richness for a first 
approximation to the problem 2 ' 16 . The binary values +1 and correspond 
to an expressed or not expressed gene, respectively. Boolean functions fi and 

the K- connection sets (9) represent the epistatic interactions among the 
genes 16 . The phenotype and/or its metabolic regulation is represented by 
the different attractors, that the dynamics of the Afif-Kauffman network gen- 
erates, each of them, playing the role of an alternative cell in the organism. 
So, the different states in the attractor represent the metabolic process 1 > 2 ' 15 . 

Changes in the fi and the mimic random biological mutations, and so, 
•d (N, K) represents the average number of genotypes giving rise to the same 
metabolic process in the phenotype by means of mutations. From the asymp- 
totic formula (26) for $(N,K), it follows that, in order that ^ represents 
a robust genotype-phenotype map; it must be a many-to-one map, which 
implies that K < K c — AK C . The number of genes that living organisms 
have, ranges from 6 x 10 3 in yeast, to less than 4 x 10 4 for the H. sapiens 14 . 
Substitution of these figures in (23) and (25) shows that, in both cases: for 
AiT-Kauffman networks, to exhibit genetic robustness; it must happen that 
K < 3 for the average number of epistatic interactions. 
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Our results are in well concordance with the fact that, the existence 
of an ordered phase representing cycles on cells, requires the emergence of 
attractors whose length grows not faster than a power of N; otherwise the 
cycle length will be too long to represent a metabolic process 2 . From the 
mean field analysis result (2) for K*, done by Derrida et. ai; that happens 
for K < 2 (in the case p = 1/2) 4 ' 10 . To remark also, is that since the 70's, 
the case K = 2 has been used as a model for cell differentiation 1,22 and the 
mitotic cycle 23 . 
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